A new integrable 3+1 dimensional generalization of the Burgers 

equation 

M. Rudnev, A.V. Yurov, V.A. Yurov 
February 8, 2008 

Abstract 

A new nonlinear 3+1 dimensional evolution equation admitting the Lax pair is presented. In 
the case of one spatial dimension, the equation reduces to the Burgers equation. A method of 
construction of exact solutions, based on a class of discrete symmetries of the former equation is 
developed. These symmetries reduce to the Cole-Hopf transformation in one-dimensional limit. 
Some exact solutions are analyzed, in the physical context of spatial dissipative structures and 
shock wave dressing. 
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1. Introduction 

The vast majority of known completely integrable nonlinear evolution PDEs is 1+1 dimensional 
[1], Uni. Dealing with more than a single spatial dimension, one faces fundamental algebraic and 
geometric obstructions. This fact accounts for the scarcity of integrable d+1 dimensional systems, 
for d >2, available today. 

A notable exception is the Burgers equation, which in 1+1 dimensions is 

Ut + UUx - VUxx = 0. (1) 

It generalizes to d > 2 spatial dimensions as follows 

+ w • Vit - z/V^tt = 0, M = -V$, V = (9^1, . . . ,9^J. (2) 

from the physical point of view, the equation describes the balance between nonlinearity and dis- 
sipation, P], rather than dispersion, which is more characteristic of integrable evolution equations. 
It comes as a model in various problems of continuous media dynamics, condensed matter physics 
cosmology, etc., see for instance [SI, [TT| . 

Mathematically, the Burgers equation is special, as it can be fully linearized via the Cole-Hopf 
(IS], |Z1) substitution 

It = -2i/Vloge, (3) 
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which reduces it to the heat equation for 0, with diffusivity i^. For the state of the art on the 
pure mathematical side of the Burgers equation, see for instance il2j and the references contained 
therein. 

The principal content of this note is as follows. Apart from the natural integrable higher dimen- 
sional generalization ((2j of the Burgers equation the latter can be used as the foundation to 
construct higher dimensional integrable PDEs, are not fully linearizable, yet their solutions can be 
found via Lax pairs. In the recent work the 2+1 dimensional BLP (Boiti, Leon and Pempinelli) 
system was studied and shown to reduce to an integrable two dimensional generalization of the Burg- 
ers equation Here we present a 3+1 dimensional nonlinear equation, which contains dissipative 
additives and has the following properties: 

1. It is a scalar second order evolution PDE with quadratic nonlinearity. 

2. In one dimensional limits this equation reduces to the Burgers equation. However, unlike the 
d+1 dimensional Burgers equation, our equation is non-isotropic, nor it is linearizable via the 
Cole-Hopf substitution. 

3. It allows for an explicit Lax pair. 

4. It possesses a class of Darboux-transformation-like discrete symmetries, and to take advantage 
of these symmetries one has to solve the Lax pair equations. The symmetries generate a 
rich spectrum of exact solutions of the equation. In particular they enable one to fulfill a 
3+1 dimensional dressing of particular solutions of the 1+1 dimensional Burgers equation. 
Conversely, in one dimensional limits the symmetries reduce to the Cole-Hopf transformation. 

2. General result 

Consider the following equation: 

K[u] =ut + ai {uj - u^x) + a2 {u-^^ - u^^) + 6i {u^Uy - u^y) + 62 {u^Uz - u^z) - pUx- 

— ^JiUy — XUz = 0, 

where u = u{t,x,y, z,t), and all other parameters are constants. 

In one dimensional limits equation Q reduces to the dissipative Burgers equation. Indeed, if 
u = u{x,t), defining 

^{x,t) = Ux{x,t), (5) 

we obtain 

6 - pCx - ai^xx + "^aiCCx = 0. (6) 
The latter equation boils down to the Burgers equation after the change t ^ t', such that 

df = dt- pdx, 

or simply letting p = 0. 

In the same fashion, if n = u{z,t) and r]{z,t) = Uz{z,t), one has 

Vt - ^Vz - a2r]zz + '2a2'nT]z = 0, (7) 

the analog of ©. 
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Finally, the reduction u = u{y,t) results in a linear equation 

Ut — ^lUy = 0. 

In view of the above, equation (jlj can be viewed as a special non-isotropic three dimensional 
extension of the Burgers equation. To emphasize this, let w = Ux, consider fi = p = X = and 
rewrite Q) as follows: 

wt + 2aiwwx + biwwy + b2WWz - aiw^x - a2Wzz - hwxy - h2Wxz - l^Wy+ , . 

+hiUyWx + b2UzWx + 2a2UzWz = 0. 

Our main result is the following theorem 

Theorem 1 Let u{x,y,z,t) be a particular solution of equation ^ and ip = 'ilj{x,y, z,t) be a 
solution of the following linear equation: 



ipt = aiipxx + a2ipzz + biipxy + b2iJxz + {p- 2aiUx - b2Uz - biUy) ipx + {p - biUx) 

+ (A - 2a2Uz - b2Ux) ipz = A[u] ip. 



(9) 



Then any Ukim = Ukimix, y, z, t), defined by the formula 

Ukim = u- log {{dx - Ux)^ {dy- Uyf {dz-UzT''^^ , {k,l,m)GZl_ (10) 

is also a solution of equation 
Theorem ^ rests on the following fact. 

Proposition 2 Equation ^ admits the following Lax pair: tpt = A{u)ijj, cf. and 



(11) 



i^xyz = Uzlpxy + Uy^^xz + Uxi^yz + {Uyz - UyUz) Ipx + {Uxz - UxUz) Ipy + {Uxy - UxUy) ■lpz + 
i'^xyz l^yz^^x '^yx'^z l^xz^^y ~l~ Ux^^yl^z) V') 

Verification of this proposition a direct calculation. Further in the note we shall refer to equations 
and (jni) as the Lax, or LA-pair for equation (jlj), and u as a potential. Observe that the spectral 
equation 1)11(1 of the Lax pair can be rewritten in a more compact form: 

{dx - Ux) {dy - Uy) {dz -Uz)'4> = -L 1 [tt] -L2 [u] ^3 M = L[u]'4) = 0. (12) 

Also observe that if we redefine the operator A\u\ A'[u] = A[u] + then the compatibility 

condition of the Lax pair equations (|1H) and © will be reduced to a an identity. Namely the 
operators L[u] and B'[u] = dt — A'[u\ will commute: [-L[n], B'[n]] = 0. 

Theorem n implies the following corollary, which follows after successive iteration of (|10() . 

Corollary 3 // {4>i}, i = 1, ••, N is a set of particular solutions of the Lax pair |3), given the 

potential u, satisfying equation new solutions of ^ are generated by the following rule: 

u = u -log (^L\^[u]L^i[u]L'^'[u]^'}j, {ki,li,m,) eZl,yi = l,...,N. (13) 

As we have indicated earlier, the Burgers equation © results from a 1+1 dimensional reduction 
of equation @. Conversely, formulae (jTU)) . (|T^ yield a bona fide generalization of the Cole-Hopf 
substitution Q. Indeed, setting n = 0, A; = / = m = Oin (|T?H) . after differentiation in x, we obtain 
the one dimensional Cole-Hopf transformation, cf. ©• Clearly, the same can be said about the z 
variable reduction as well. 
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Proof of Theorem [T] 

The theorem wih fohow from the fohowing lemma. 

Lemma 4 Let ip be a solution of the Lax pair equaitons Ul\) . ^ with the potential u, which is a 
solution of equation 0). Then the function 

ijkim = L\[uWML'^Hi^. {k,l,m) £Zl (14) 
also satisfies Ul\) . with the same potential u. 

To prove the lemma, observe the validity of the following commutator relations, for i,j = 1,2,3: 

[U[ulL[u]] = [Li[u],B[u]] = [L^[u],Lj[u]] =0. 

Lemma 0] then follows by induction. □ 

Now, to prove Theorem ^ let us introduced three intertwining operators 

Di = fidi-gi, i = 1,2,3, (15) 

with the quantities fi,gi to be found (naturally, i9i,2,3 = dx,y,z, respectively), such that the operators 
Di have the following property: for some Ui = Ui{x,y, z,t), 

L{ui)D, = DiL[u], B{ui)Di = DiB[u]. (16) 

The commutation relations H16() determine the maps u ^ Ui, which come from substitution of H15|) 
into (|l(i|l . The explicit form of the operators Di is found as follows. 

Substituting (|15() into ()16|) and equating the components at the same partial derivatives results 
in a system of nonlinear equations (which is not quoted because of its bulk) whence it follows: 

Di = e''' {Li[u]- a) , Ui = u = u - V, (17) 

where Cj - are constants, which will be further assigned zero values. The quantity v = v{x, y, z, t) is 
a solution of the following nonlinear equation: 

Vt = ai [v^x + Vx"^) + 02 (^22 + Vz'^) + hi {Vxy + VxVy) + 62 {Vxz + VxVz) + 

+ {p- 2aiUx - h2Uz - huy) Vx + {fJ, - hux) Vy + {X - 2a2Uz - b2Ux) Vz 

Therefore (fT7]) or explicitly (|T8|) indicate that for n = 0, the function "—v" satisfies equation Q). 

Then automatically the quantity = e"" will satisfy the Lax pair equations ()11() and @. In 
fact, the L-equation, cf. (fTT|) . is satisfied as the identity. The A-equation Q however is satisfied 
only if ?x is a solution of (@J). On the other hand, by Lemma^ the functions of the form ipkim defined 
via relation H14|) are also solutions of (|11() and ©, with the same potential u. Rewriting them as 
'ipkim = exp(t'fcjm) and substituting into © one verifies that the quantities Vkim are indeed solutions 
of equation ()18|) . Theorem ^ and formula (|1U() now follow from the second relation from H17|) . □ 

Remark. Formula (|10)) has a countenance similar to the Darboux transformation, which is a stan- 
dard tool for construction of exact solutions of nonlinear PDEs (usually 1+1, more rarely 2+1 
dimensional) which admit Lax pairs, see e.g. IR for the general theory, applications and references. 
However ()1U() does not represent a bona fide Darboux transform for two following reasons. 
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1. Darboux transforms, representing discrete symmetries of a Lax pair, possess a non-trivial 
kernel on the solution space of the pair. In other words, there always exists some Lax pair 
solution which zeroes the transform. This is the property which enabled one Crum, jH], to 
write down the determinant formulae for successive Darboux transforms. Transformation H1U|) 
however does not have this property. 

It is known that in addition to the Darboux transform, fairly rich spectral problems, such as 
the Zakharov-Shabat problem for the Nonlinear Schrodinger equation or its two-dimensional 
generalization for the Davey-Stewartson equations, admit another discrete symmetry, namely 
the Schlesinger transform, |H]. The difference between and the latter transformations lies 
in the fact that for the Schlesinger transform, the potential transformation rules can be locally 
defined without using the Lax pair solution ^, while (|lflj) certainly does so. This feature is 
shared by ()1U() and the standard Darboux transformation. 

2. To construct exact solutions of nonlinear PDEs via the Darboux transform, one has to take 
advantage of the solution of the full Lax pair as a system of equations. In order to get (|Tn|) 
however, we have used the solution of the A-equation Q only. The L-equation ((TT|) has been 
used only as a tool to prove Theorem ^ To this effect, transformation combines essential 
features of the Darboux and Cole-Hopf transformations. 



4. Some exact solutions 

Let us use the above formalism in order to construct some exact solutions of equation (jH} (which 
is the equation for Ux, where n is a solution of equation (jlj) with fj. = X = p = 0). We consider 
equation ©, because it appears to be a closer relation of the Burgers equation and is likely to be 
interesting from the physics point of view. 

Example 4.1. As the first example let us consider dressing on the vacuum background n = 0. 
In this case the function 

ipix, y, z, t) = a^x^ + h^y"^ + (?z^ -h 2 (oia^ + a2(?) t ■\- , (19) 

where a, 6, c, s are some real constants, is clearly a solution of the Lax pair equations (|11() and ® . 
Substituting (tTTHl in (fT(l we derive the Ukim- After differentiating it with respect to x and choosing 
k = I = m = 0, we obtain a solution w of equation © as follows: 

w{x, y, z, t) = - ^2^2 ^ ^2y2 + ^2^2 + 2 {aia^ + a2C^) t + s^' ^ ' 

Physically, solution (|2()j) describes a rationally localized impulse, vanishing as t — > -j-oo. To ensure 
that H2U() is non-singular for t > 0, one should impose the inequality 02 > —a^axjc? on the coeffi- 
cients. Moreover, if 02 = —a^ai/c?, the solution in question is rationally localized and stationary. 

Remark. The fact that there exists a localized stationary solution in an equation containing 
dissipative terms may appear somewhat paradoxical from the point of view of physics. However 
solution ()2U() is stationary only if aia2 < 0. One can see that these constants appear in the dissipative 
terms of equation Q, ai characterizing the dissipation along the x and 02 along the z axes. The 
fact that ai and 02 should have different signs implies that dissipation in the direction of one axis 
is compensated by instability in the direction of the other. The balance of these two effects results 
in the stationary solution, which can be regarded as a three-dimensional dissipative structure. A 
similar situation occurs with two-dimensional stationary solutions of the BLP equation, cf. |13j . 
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Example 4.2. Another solution of (jlH) . @ when n = is 

tpix, y, z, t) = cie"(""i+^''i)* cosh(ax + l3y) + c2e('^i'^'+'^262+afef>2)t cosh(ax + bz) + cge'^^c^i cosh(cz), 

(21) 

where a, /3, a, b, ci , C2, C3 are some real constants. Choosing them such that 

(3=-^, f^^zh±ME^ (22) 
bi 2a2 

and using and (|22|) in the same way as in Example 4.1 above (|2U|1 . we obtain another solution: 

aci sinh(ax + /3y) + ac2 sinh(aj; + bz) 

w{x, y, z, t) = — , —2^ 

ci cosh(ax + py) + C2 cosh(ax + bz) + cae'^^c t cosh(czj 

Example 4.3. Developing the analog with jl.Sj further, one can develop a procedure of con- 
struction of exact solutions of the three-dimensional equation Q which are based initially on the 
solutions of the 1-1-1 dimensional Burgers equation. Consider equation © for the unknown ^. 
Suppose A = ;U = p = 0, let oi = in @. 

Clearly the quantity U{x,t) = 2v^[x,t) solves the one dimensional Burgers equation 

Ut + UUa, - uU^^ = 0. (24) 

As a starting point let us take a shock wave solution of (|24|1 . e.g. 

V — va a 

^ = ^^ = ^ + TT^^' 

where a and v are constants. Seek a solution of © as a superposition 

N 

^ = ^^fc(ry)e/3'=y+T'=^, (26) 

k=l 

where rj = x — vt, while the 2A^ quantities /3fc and 7^ can in general be functions of ry. For simplicity 
however let us render them constants to be determined. 

Substitution of into © yields N linear equations for Ak{ri): 

uA + {v + a- 2vC) A + (027^ -ai)A = 0. (27) 

In equation (P7|) the indices for the quantities Ak, (3k, 7fc and crfc = bi(3k + b2^k have been omitted, 
while the quantity came from (|25() . dot standing for differentiation with respect to r/. 

Equation H27|l can be simplified further. To do so, let us introduce a new independent variable: 

^/ N /- on o ^ ^ V — va 

q = i{x,t)-io, e'^^ = --l, ^0 = ^^ . 

q 2v 

In terms of q equation \21\ becomes 

V {q^ - aqf^'{q) + a {q^ - aq) ^ {q) + {5- aq) A{q) = 0, (28) 
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where 6 = 027^^ — ct^q, while prime stands for differentiation in q. The latter equation can be 
simplified further via a substitution 

Aiq) = Wiq) M- , (29) 



reducing to the following equation for W{q): 

W"{q) _ 2aav - 4.5v + a"^ 
W{q) ~ Au'^iq-a^q^ ' 

In a particular case of the dependence 



(30) 



h = r ; 

h 

one can solve ^ explicitly as 

N / „ \ '^k/i'iya) 



k=l ^ 



q — a 



(31) 



where and are arbitrary constants. Substitution of (|5T|) into ((TU)) yields exact solutions of 
equations @, ©• The described procedure can be regarded as the shock wave dressing. 
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